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ABSTRACT 

Fermi Gamma ray Space Telescope observations of the flat spectrum radio 
quasar 3C 454.3 show a spectral-index change AF = 1.2 ± 0.3 at break energy 
Ebr ~ 2.4 ± 0.3 GeV. Such a sharp break is inconsistent with a cooling electron 
distribution and is poorly flt with a synchrotron self-Compton model. We show 
that a combination of two components, namely the Compton-scattered disk and 
broad-line region (BLR) radiations, explains this spectral break and gives a good 
flt to the quasi-simultaneous radio, optical/UV, X-ray, and 7-ray spectral energy 
distribution observed in 2008 August. A sharp break can be produced indepen- 
dent of the emitting region's distance from the central black hole if the BLR has 
a gradient in density oc R~^, consistent with a wind model for the BLR. 

Subject headings: gamma rays: galaxies — Quasar: 3C 454.3 — radiation mech- 
anisms: non-thermal 



1. Introduction 

The flat spectrum radio quasar (FSRQ) 3C 454.3 (PKS 2251 + 158) is representative of 
blazars with large, > 10^^ erg s~^ apparent 7-ray luminosities, and is one of t he brightest 



sources seen with the recently launched Fermi Gamma-Ray Space Telescope (iTosti et al. 



20081 : lEscande fc Tanakal 120091 ). It has also been detected with AGILE up to ^ 3 GeV , 



including a giant outburst in 2007 prior to the Fermi launch ( iVercellone et al.ll2009l . |2010| ). 
Based on > 100 MeV data taken during the checkout and early science phase of the mission 
(2008 July 7- October 6), this source was found to be variable on timescales down to ti, ~ 



-2- 



few days (lAbdo et al.l |2009| ). A well-sampled data set using the ground-based SMARTS 
(Small and Moderate Aperture Research Telescope System) in the B, V, R, J, and K bands, 
and the Swift X-ray telescope (XRT) and ultraviolet-optical telescope (UVOT) was col- 
l ected simultaneously with Fermi observations in the period 2008 August through December 
( iBonning et al.l 120091 ) . These multiwavelength observations revealed that the infrared, op- 
tical, and 7-rays were have similar variability patterns. The 7 rays display similar relative 
flux changes compared to the longer optical w 'avelengths, but small er relative flux changes as 
the short wavelength optical emission, which iBonning et al.l (120091 ) interpreted as being due 
to an underlying accretion disk. The Swift X- ray d ata did not strongly correlate with the 
other wavelengths, though iJorstad et al.l ( 120101 ) and iBonnoli et al.l (120 lOl ) flnd a correlation 
at different epochs. 

The Fermi observations of 3C 454.3 revealed a broken power-law spectrum, with F = 
2.3 ± 0.1 (where the photon flux $ oc E~^) and F = 3.5 ± 0.25, below and above a break 
energy of 2.4 ± 0.3 GeV, respectively JAbdo et al.lbood l This spectrum is not consistent 
with Compton scattering of a lower-energy photon source by jet electrons in the fast- or slow- 
cooling reg imes, the sta i idard interpretation for the high energy spectrum in leptonic models 
of blazars. lAbdo et al.l (120091 ) also consider it unlikely that the spectral break is a result of 
77 absorption with a local radiation fleld or the extragalactic background light (EBL), the 
former because this would require a lower jet bulk Lor entz factor than tha t inferred from 
superluminal radio observations of 3C 454.3 {Ti,uik ~ 15; IJorstad et al.ll2005l ). and the latter 
because the universe is transparent to 40 GeV photons a t 2: 5, 1-0 for all EBL models (e.g., 
Finke et al.lboiol : Istecker et allbood : Icilmore et alll2009[ ). lAbdo et all J2009[ ) thus consider 
the most likely explanation to be that the spectral break is caused by an intrinsic break in 
the electron distribution. However, as explained below, this is unlikely, given the broadband 
spectral energy distribution of this source. 

In this letter, we show that a combination of two Compton scattering components can 
explain the break in the LAT spectrum of 3C 454.3. In particular, we consider a combination 
of components from Compton-scattered accretion disk and broad-line region (BLR) radia- 
tion, both of which should be strong seed photon sources in FSRQs. The implications of this 
model for the BLR and for the variability of the emitted radiation are explored. Finally, we 
conclude with a brief discussion, with implications for breaks in other sources. 

We use parameters Hq = 71 km s~^ Mpc"'^, Qm = 0.27, and = 0.73 so that 3C 454.3, 
with a redshift of z = 0.859, has a luminosity distance of di = 5.5 Gpc. 
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2. The Spectral Energy Distribution of 3C 454.3 



The contemporaneous SED data from 3C 454.3 are given in lAbdo et al.l (l2009f ) and 
plotted in Fig. [1] The radio data are likely from a larger region of the jet than the rest of the 
SED, and can be treated as upper hmits for the purposes of model fits. We assume that the 
electron distribution which produces this radiation has the form of a broken power-law given 
by rie oc 7'"^' for i^^^ < 7' < 7^^^ and Ue oc 7'"^'^ for 7^^^ < 7' < i^^^. The Swift UVOT data 
give a power-law spectrum with T^pt ~ 2.9, presumably from electron synchrotron emission. 
The X-rays from the Swift XRT give Tx ~ 1-4, which should be the index of the synchrotron 
emission below the break, assuming that the X-rays are from Compton-scattering with the 
same electrons as are causing the lower-energy synchrotron emission. Thus the Swift data 
imply that, if the electron distribution can indeed be represented by a double power-law, 
Pi ~ 1.8, and p2 ~ 4.8, since p = 2r — 1. This spectrum indicates a fast cooling regime, 
where one expects pi ^ 2 and p2 = Q' + 1, where q is the injection index. 

If the LAT spectrum was the result of Thomson scattering of a single source of photons, 
whether that source is the disk, BLR, or dust torus, the spectral indices would be about the 
same as for the synchrotron, namely, F 1.4 and F 2.9 well below and above the break, 
respectively. Instead, F ^ 2.3 and F ^ 3.5 are observed, respectively, although it is possible 
that Klein-Nishina (KN) effects could soften the higher energy index. The lower energy LAT 
index furthermore does not necessarily correspond to the asymptotic value at ~ E^j.. 



Abdo et al.l (|2009[ ) suggested that the break could be caused by the high energy cutoff 
in the electron distribution, 7maa;- However, if this is the case, below the 7-ray break, one 
would expect the same F as in the optical {Topt = 2.9), because of their correlated variability, 
instead of the F = 2.3 which is observed. The 7-ray spectral index could be softer than Topt, 
due to KN effects, but it is unlikely to be harder. 

It is possible that the X-rays through Fermi 7-rays are all created by a single component, 
likely the SSC component, since this is generally broader than external Compton components. 
In this model, the softer 7-rays are created by the portion of the SSC component at and 
below the break, as it gradually softens. This SSC-only model gives an adequate fit to 
the broadband SED, although it gives a poor fit to the LA T data (Fig , [ll), a nd does not 



produce a sharp break in the LAT 7-ray spectrum, although I Abdo et al.l ( 120091 ) found that 
the LAT spectrum could just as easily be fit by a smooth transition such as, for example, a 
power-law w i th an exponential cutoff. See Table 1 for the spectral parameters for this fit and 



Finke et al.l ( l2008l ) for details on this model. The SSC-only model in Fig. [T] does not enter 
the KN regime until around ^ 3 x 10^^ Hz, which is considerably above the Fermi-LAT 
data, and these 7-rays are thus produced entirely in the Thomson regime. However, it is 
possible that local 77 absorption could soften the higher energy portion of the spectrum and 
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produce a better fit. 

But the SSC model has other problems. It ignores the BLR for this source, which must 
be present due to the strong broad emission lines seen in this object's optical spectrum; it 
results in an extremely low magnetic field, B, which is considerably lower than the equipar- 
tition value, Beq, and thus the particle jet power, Pj^par is much greater than the magnetic 
field jet power, Pj ^ and close to the Eddington limit (see Table 1); and it may not be con- 
sistent with the variability of th is source (see § below). The constraints on the Doppler 



factor from iJorstad et al.l ( l2005l ) and the constraints on the emitting region size from the 
variability timescale imply B = 0.032 ± 0.005 G for the SSC model. We propose that the 
LAT spectrum of 3C 454.3 is instead the result of a combination of two Compton-scattered 
components, for example, Compton scattering of the disk and BLR radiation. 



3. BLR Radiation Energy Density 

Before modeling the SED of 3C 454.3, we derive approximate formulae for the energy 
density of the BLR at the location of the jet. A regime where the energy densities of the 
Compton-scattered disk and BLR components are approximately equal can explain the LAT 
spectrum of 3C 454.3. 

We assume that the BLR is spherically symmetric and extends from inner radius Ri to 
outer radius Ro, and the electron number density of the BLR is a function of distance from 
the central source, nBLR^R) = no{R/RiY , so that the Thomson depth of the BLR is given 
by tblr = c^T Jr" dR nBLR{R) giving, when Ri < Ro, 




c<-i 



The BLR Thomson-scatters photons from the accretion disk, assumed to be a point 
source located at origin with luminosity Ld- If the distance from the black hole to the 
7-ray emitting region, r < Ri, the energy density from BLR radiation is approximately 
independent of r, so 



„ , LrfCXT dR 
UBLR[r < Ri) = j ublrKR) 

C<-i 



1 c±i 

T I "^^"^ ' . 

^dTBLR J 1 f R,y (^+1 



AUC 1 RiRo \Ro I 

1 C±i 



-1<C<1 ' (2) 
C>1 




V [Hz] 

Fig. 1. — The SED of 3C 454.3. The Effelsberg/IRAM radio data are given by the diamonds, 
the Swift-UVOT and XRT data are given by the squares and triangles, respectively, and the 
Fermi-LAT data by the filled circles. The model components for the multi-component model 
are shown as dashed curves, where the purple dashed curve is the Shakura-Sunyaev disk, 
the red is the synchrotron, the green is the SSC, the blue is the Compton-scattered disk, 
and the brown is the Compton-scattered BLR component assuming tblr = 0.01. The thick 
dotted curve gives the SSC fit, and the thick solid curve shows the multi-component model. 
The inset shows the Fermi-LAT spectrum in more detail. The arrows indicate the frequency 
of the spectral break. 
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Table 1. Model Parameters^ . 



Symbol 


Multi-Component model 


SSC model 


z 


0.859 


0.859 


^bulk 


15 


15 


Sd 


30 


15 


B 


0.4 G 


0.032 G 


tv 


1.7 X 10^ s 


1.7 X lO'"^ s 


PI 


2.0 


1.8 


P2 


4.4 


4.8 


7' ■ 


3 X 10^ 


10 




2 X 10* 


2 X 10^ 




1.1 X 10^ 


10* 


A/9 


2.0 






0.04 




V 


1/12 




r 


1.5 X lO^Rg 




TBLR 


0.01 




R^ 


5.0 X lO^flg 




Ro 


5.0 X lO^Rg 




C 


-2 




B/ Beq 


0.6 


0.06 




1.8 X 10*^ erg s"! 


10*^ erg 


p 

J,par 


2.7 X lO*** erg s"! 


2.3 X 10*^ erg s"! 



^See text and lDermer et aL for details on parameters 

and calculations of jet powers. 
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usmg eqn. 



Dermer et al. 



j)^ If Ri < r < Ro, then, assuming the logarithm term in eqn. (93) of 

J2009h is ^ 1, 



UBLR[r) 



Ancr 



Ro 



dR 
~R 



nBLR{R) 



47rc 



1 / J_ 

r-Ri \ Ri 
1 

rRo \ Ro 



C+1 

c 

' C±i 

c 



1 



C+1 



RiRo y Ro j C, — \- 

1 c±i 



(see also lDonea &: Protheroell2003l ). 



c<-i 

-i<C<o 
< c < 1 



(3) 



4. Compton Scattering Disk and BLR Radiation 



The external Compton-scattered BLR vF^, flux in the Thomson regime can be approx- 
imated as 



ECBLR 



CaTUBLRj'^Neij 



(4) 



where 7' = ^/ e{l + z)/ [5^6^(8. 4rg)] ( IDermer fc Schlickeiserl Il994l ) a nd e is the observe d 
photon energy. For 3C 454.3, T^y^^ik ~ 15, the Doppler factor 5^ ~ 25 (jjorstad et al.ll2005l ). 
Gu et al.l J200lh found Mg = MB/f/(109MQ) ^ 4.4 where M^ h is the black hole mass. A 
more recent analysis of the black hole mass gives Mg = 1 — 3 (IBonnoli et al.l I2OIOI ) , and we 
assume for the remainder of this work that Mg = 2. If the 7-ray emission region is located 



at r < Tl^if^Tg ^ 



10^ pc, where r„ = GMbhc 



1.5 X 10 Mg cm, then the external 



Compton-scattered disk emission is in the near field regime. In this case the received vFy 
flux can be approximated by 



ECD 



5^ caT 



(5) 



where f = 2^(1 + z)e/ [5led{\/?>r)\ (IDermer fc Schhckeisei]l2002h . Here 

lEdd 



ed{R) = 1.5 X 10" 



n 1/4 \ -3/4 



MgT] 



is the mean photon energy from a IShakura fc SunyaevI (Il973l ) accretion disk at R, Isdd 
Ld/LEdd, LEdd = 1-3 X lO'^^Mg, and r] is the accretion efficiency. 
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A solution where the Compton-scattered disk and BLR emission are approximately 
equal can explain the Fermi-LAT emission. Also, a robust solution is obtained by noting 
that both f^'^^^R and f^^^ are oc r"^ if C = -2 and Rj < r < Rn- A BLR gradient o f 



C = —2 is consistent with a wind model (IMurray fc Chianglll995l : IChiang fc Murraylll996l ). 



Setting eqns. @ and ([5]) equal to each other, and, ignoring the fact that the electrons 
generally scatter photons of different energies, one obtains the condition tblr = (Ri/f^g)'^ 
for emission regions formed in the BLR. In order to have comparable flux from the two 
components, this relation shows that the total Thomson depth through the wind should 
increase in inverse proportion to the inner radius of the BLR. 



An accurate fit to the multiwavelength SED using the mult i- component model of lDermer et al 



( 120091 ). using model parameters from Table 1 with t^lr = 0.01, is seen in Fig. [TJ In this 
model, the portion of the Fermi spectrum below the break is from Compton scattering of a 
combination of direct accretion disk and BLR-reprocessed radiation, while above the break 
it is almost entirely from Compton scattering of BLR-reprocessed disk photons. This model 
uses the full Compton cross section for all Compton calculations, and the full Shakura- 
Sunyaev disk spectrum for both direct and BLR-reprocessed Compton-scattered emission. 
For the Compton-scattered BLR calculation, the Shakura-Sunyaev disk spectrum is assumed 
to originate from the black hole itself, whereas in the direct disk-scattered calculation, geo- 
metric effects are taken into account. 

For KN effects to become important in Compton scattering the BLR radiation, 7' > 
{4Ti,uik^d)^^ where is the typical dimensionless energy of photons from the disk in the 
disk's stationary frame. For our model, = 2.2 x 10~^, and Tbuik = 15, so that 7^^^ = 760. 
Electrons with this energy scatter photons primarily to energies of exN ~ '^D^<i7^Ar/(l + -2) = 
6.2 X 10^ which corresponds to frequencies oi u = 7.6 x 10^^ Hz. Thus above the break the 
Compton scattering of the BLR is almost entirely in the KN regime. Since the KN cross 
section goes as approximately ckn e~^, ff^^ should go as F Topt + 1 = 3.9 above 
the break in the 7-ray spectrum. This is quite close to the F = 3.5 which is observed by 
Fermi. Although the BLR photons are Thomson-scattered disk photons with essentially 
the same energy as the direct disk photons, geometric effects from Doppler boosting and 
aberration make Compton-scattered BLR radiation peak at higher photon energies than the 
disk photons, which tend to come from behind. 

The threshhold for 77 absorption from the BLR-reprocessed disk radiation implies that 
7-rays with energies e > 2[erf(8.4rg)(l -|- z)]~^ will be absorbed. For our model this turns out 
to be e ~ 4.9 x 10^ ot u ^ 6.0 x 10^*^ Hz, which is above the highest energy photon observed 
with the Fermi-LAT from this source. 

The sum of the particle and magnetic- field power Pj^par + Pj,B is well below LEdd = 
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2.6 X 10^''' erg s ^ for Mg = 2.0 (the photon power is ~ 8 x 10^^ erg s ^). The model is close 
to equipartition between the total particle and magnetic-field energy densities, assuming a 

10R>. Our 



J, par j,f,- 



factor 10 X more energy in total nonthermal particles than electrons, i.e., Pj^ 
multi-component fit has lower jet power than the fit to the same SED by 
(120091 ) ■ due primarily to their assumption of a larger proton power, with Pj^par > 200 x Pj ^. 



Ghisellini et al. 



5. Variability 



The 7-ray spe ctral break observed with Fermi was constant over several month timescales 
( lAbdo et al.ll2009l ). As the 7-ray-emitting region moves outward, the spectral break should 
remain at approximately the same energy, since for ( = —2, the BLR- and disk-scattered 
fiux should remain approximately equal. The overall 7-ray fiux will, of course, change over 
time, but as long as i?i < r <^ T^r^, and Ri or tblr does not change significantly, the 7-ray 
spectral break should be present. 



In iBonning et al.l (|2009[ ). optical observations are seen to vary with changes on the 
same magnitude as the 7-rays. If the 7-rays were produced by SSC emission, "quadratic 
variability", i.e., f^'^'" oc (Z^^")^, might be expected. The absence of such variability is in 
accord with our spectral fitting results means that this model for 3C 454.3 is disfavored. 

The X-ray light curve of 3C 454.3 shows little variability duri ng this time, and is not 
correlated with the IR, optical, UV, or 7-rays (jBonning et al.l 120091 ). This can be explained 
by the longer cooling timescale (see Fig. [2]). These timescale calculations use monochromatic 
approximations for the Compton-scattered disk and BLR components, but the full Compton 
cross section for all components. The electrons which make the X-rays by SSC have lower 
energies than those that make the optical and 7-rays, and thus cool more slowly than the 
observed variability timescale. 

The lack of correlated variability between X-rays and optical and 7-rays might indicate 
that the X-rays originate from another part of the jet than the LAT 7-rays. However, 
in this case, the emission would still likely be in the fast cooling regime, since the 7-ray 
fiux is significantly greater than the optical and radio fiuxes, and as long as one is in the 
fast cooling regime, Compton scattering of a single component cannot explain the LAT 
spectrum. Furthermore, during the 2007 November fiares of 3C 454.3, detected by Swift 
UVOT, XRT and BAT, INTEGRAL, and AGILE, the X-rays were we akly variable, compare d 
with other wavebands, but did seem to show correlated variability ( IVercellone et al.ll2009l ). 
which indicates they originate from the same emission region. One final possibility is that a 
hadronic model may be possible to explain the LAT 7-ray spectrum (e.g.. lMucke &: Protheroe 
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Fig. 2. — Cooling timescales in the observer's frame. The electrons with 7 which pro- 
duce the majority of the radiation in certain wavebands are marked. The electrons which 
produce X-rays from SSC emission have total cooling timescales longer than the variability 
timescale, explaining the lack of X-ray variability. The optical synchrotron radiation is made 
by electrons with 7 ^ 6-8x10'^, which would therefore correlate with the 7 rays formed by 
Compton-scattered disk radiation. 
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200ll ). however, such models may have difficulties in explaining the correlated variability seen 



in this source, since protons would evolve on longer timescales than electrons. 



Discussion 



Assuming the dual-component Compton-scattering scenario presented in this paper is 
correct, what can it tell us about the location of the 7-ray emission region? It implies that 
this region is < r ^ T^r^, so that, using the parameters from the modeling of 3C 454.3, 
the location of the jet falls over a large range from 0.1 pc < r ^ 100 pc. However, the 
external energy density from disk and BLR radiations is oc r~^, a steep decline, so that 
the blob must be within 0.1 pc for a significant scattering component. Recent modeling 

201c) also finds that the 7-ray emitting region is within the 



of 3C 454.3 flBonnoh et al 



BLR. By contrast. 



Sikora et al. 



(120081 ) suggested that the 7-rays observed by EGRET and 
AGILE were caused by Compton-scattering of radiation from hot dust by a blob located 
~ 10 pc from the central black hole, though detailed modeling is still needed to determine 
if this model can explain the sharp spectral break in the Fermi observations. The blazar 3C 
279 also shows a similar spectral break during the onset of an extended episo de of optical 



polariz ation, while the jet emission region could still be within the BLR region ( lAbdo et al. 
2010ah . 



Whether this model can predict similar spectral brea ks in other FSRQs and low-synchrotron- 
peaked blazars, such as AO 0235+164 or PKS 1502+106 JAbdo et al.lboiObh . depends mainly 
on the properties of the BLR in these sources, in particular, whether tblr ~ ' '^nlRj- spe- 
cific blazars, especially at higher z, 77-absorption by scattered disk radiation (jReimeiil2007l ) 
could still be effective. Detailed modeling of simultaneous SED data will test these scenarios, 
as it does for the c ase of 3C 454.3. If our model is correct, it provides evidence for a wind 
model of the BLR jMurrav fc ChianeJll995l : Ichiang fc Murravlll996l : IeIvIsI boooh . 
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